Gluon and quark jet multiplicities at N 3 LO+NNLL 
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We present a new approach to consider and include both the perturbative and the non-perturbative 
contributions to the multiplicities of gluon and quark jets. Thanks to this new method, we have 
included for the first time new contributions to these quantities obtaining next-to-next-to-leading- 
logarithmic (NNLL) resummed formulae. Our analytic expressions depend on two non-perturbative 
parameters with a clear and simple physical interpretation. A global fit of these two quantities 
shows how our results solve a longstanding discrepancy in the theoretical description of the data. 
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Collisions of particles and nuclei at high energies usu- 
ally produce many hadrons. In Quantum Chromody- 
namics (QCD) their production is due to the interac- 
tions of quarks and gluons and to test it as a theory 
of strong interactions, the transition from a description 
based in terms of quarks and gluons to the hadrons ob- 
served in experiments is always needed. The production 
of hadrons is a typical process where non-perturbative 
phenomena are involved. However, the hypothesis of Lo- 
cal Parton-Hadron Duality (LPHD) assumes that parton 
distributions are simply renormalized in the hadroniza- 
tion process without changing their shape allowing 
perturbative QCD to make predictions. The simplest ob- 
servables of this kind are gluon and quark multiplicities 
(rih)g and (n/ l ) s which represent the number of hadrons 
produced in a gluon and a quark jet respectively. In 
the framework of the generating-functional approach in 
the modified leading logarithmic approximation (MLLA) 
0, several studies of the multiplicities have been per- 
formed In such studies, the ratio r = (nh) g /{rih)s 
is at least 10% higher than the data or it has a slope too 
small. Good agreement with the data has been achieved 
in Ref. 6] where recoil effects are included. Nevertheless 
in Ref. 6] a constant offset to be fitted to the quark and 
gluon multiplicities has been introduced, while the au- 
thors of Ref. [7[ suggested that other, better motivated 
possibilities should be studied. 

In this Letter, we study such a possibility inspired by 
the new formalism that has recently been proposed in 
Ref. Thanks to very recent new results in small- x 
timelike resummation obtained in Ref. [s| , we are able to 
reach the NNLL accuracy level. A purely perturbative 
and analytic prediction has been already attempted in 
Ref. 0] up to the third order in the expansion parameter 
tJHI i.e. a s , where paradoxically the quark multiplic- 
ity and the ratio are not well described even if the be- 
havior of the perturbative expansion is very good. Our 
new resummed results that we present here are a gener- 
alization of what was obtained in Ref. 0| and represent 
also a solution to this apparent paradox. 

We consider the standard Mellin-space moments of the 



coupled gluon-singlet system whose evolution in the scale 
is governed in QCD by the DGLAP equations: 
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The timelike splitting functions Py can be computed per- 
turbatively in the strong coupling constant: 
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where u = N~ 1 with N being the usual Mellin conjugate 
variable to the fraction of longitudinal momentum x. The 
functions P^\w) with k = 0, 1,2 appearing in Eq.@ in 



the MS scheme can be found in Refs. [IM3 through 
NNLO and in Refs. H E! Gil through the NNLL. 

In general it is not possible to diagonalize Eq.([T]) be- 
cause the contributions to the splitting function matrix 
do not commute at different orders. It is, therefore, con- 
venient (see e.g. Ref. [Hj]) to introduce a new basis, called 
"plus-minus" basis where the LO splitting matrix is di- 



agonal with eigenvalues P|°| and P^'l . We define such a 
change of basis according to the following transformation 
of the gluon and the singlet fragmentation functions in 
Mellin space: 
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where 
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The general solution to Eq.([T]) can be formally written 
as 

P(/i 2 ) = jexp^ ^P(f)\ D(f4), (5) 
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where denotes the path ordering with respect to /i 2 
and 
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Now making the following ansatz: 
^ dfi 2 1 



T 



exp 



t>r, 



where 



Z _1 (M 2 )exp 



dp 2 



P ++ (u>) 
P- 



Z(Mo), (7) 
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is the all-order diagonal part of the splitting matrix in 
the "plus-minus" basis and Z is a matrix in the same 
basis with a perturbative expansion of the form: 



Z^ 2 
we obtain that 
D a (cj,fi 2 ) = D+(cj,fi 2 ) + D-(cj,fi 2 ); 



l + aV^W + Ota 2 ), (9) 
a = s,^, (10) 



where Z? + (w,/i 2 ) evolves like a "plus" component, 
D~(u>,fj, 2 ) evolves like a "minus" component, and 

^,,x 2 ) = bt{u,ti)T ± {u,^,nl)Ht{uj^ 2 ). (11) 

Here T±(u>, [i 2 , [Iq) is a renoramlization group exponent 
which is given by 



f±(u,n 2 ,nl) = exp 
with 



aAf * 2) da. 



a s (>5) 
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P±±(u),a s ) 



, (12) 



dfi 2 

Poa 2 s (ti 2 )-Pia*(p 2 ) + O(ai). (13) 



We recall that 



Po = —C A ~ 3 n / T R' 
34 20 

fh = —C 2 A -—C A n f T R -4C F n f T R , (14) 

where C A = 3, C F = 4/3 and T R =l/2 in QCD and n/ 
is the number of active flavors. In Eq. fTTj) H^(uj, /x 2 ) are 
perturbative functions containing off-diagonal terms of P 
beyond the LO and the normalization factors D^(uj, /Jq) 
satisfy the following conditions: 



D+{u,4) = — — — Dg[ui, /j,q) ; 



D s (w,Mo) 



(15) 



We note that differ from Df starting at higher orders 
(lij . In the following, we collect the resummed formulae. 
Details of the calculation will be presented elsewhere [l6[ . 

After the resummation is perfomed fo r P++ i n Eq.© 
thanks to the results obtained in Refs. (9l.ll7lll8|. we find 
for the first Mellin moment (u> = 0) at NNLL: 



P™ ll (oj = 0) = 7 o(l - A' l7o + K 2 ^ + O( 7o a )), 
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with £2 = 7r 2 /6. Now we can perform the integration in 
Ea.([I21) up to the NNLL to obtain that 
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We are now ready to define the avarage multiplicities 
in our formalism: 

(n h (Q 2 )) a = D a (0,Q 2 ) = D+(0,Q 2 ) + £-(0,Q 2 ), (23) 

with a = s,g for the quark and gluon multiplicities re- 
spectively. From Eqs. (|llll5|) we have that 



-D+(0,Q 2 ) 
and 

Ds{0,Q 2 ) «™ 



lim 



aui H+(Lj,Q 2 ) 



l-g„ff g ->,Q 2 ) 
e w Hs{uj,Q 2 ) 



r+(Q 2 ), (24) 



= r_(Q 2 ). (25) 
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Using these definitions, it is convenient to write for the 
gluon and quark multiplicities in general: 

(n h (Q 2 )) g = D+(0,Q 2 )f r + es (0,Q 2 ,Q 2 )H+(0,Q 2 ) 



D-(0,Q 2 )r^(Q 2 )f r J s (0,Q 2 ,Q 2 )H-(0,Q 2 ) 7 
fl es (0,Q 2 ,Q 2 )H+(0,Q 2 ) 



{n h {Q ))s = 



r+{Q 2 ) + ^ g 
+D-(Q, Ql)t: es {Q, Q 2 , Q 2 )H- (0, Q 2 ). (26) 

For the coefficients of the renormalization group expo- 
nents, we clearly have the following simple relations at 
the lowest order in a, 



r + {Q 2 )=C A /C F - 
Hf{Q 1 Q 2 ) = l- 



MQ 2 )=0; 
Dt(0,Ql) = Dt(0,Q 2 ),(27) 



with a — s, g. One would like to include higher-order cor- 
rections to Eq. (|27p . However, this is highly non-trivial 
because the general perturbative structures of the func- 
tions if„ and Z±zp^ a , whose knowledge is required for 
the resummation, are not known. Fortunatly, general as- 
sumptions and approximations can be made to improve 
them. Firstly, it is a well known fact that the "plus" com- 
ponents by themselves represent the dominant contribu- 
tions for both the gluon and the quark multiplicities (see 



e.g. Refs. [19|,|20|). Secondly, Eq.([25]) tells us that D~ is 
suppressed with respect to D~ , because a w ~ 1 + 0(uj). 
These two facts suggest us that to keep r_(Q 2 ) = even 
at higher orders should still represent a good approxi- 
mation. Then we notice that higher-order corrections to 
Dq(Q, Qq) and H^(0, Q 2 ) just represent a redefinition of 
Dq(0, Qq) apart from running coupling effects starting 



at order 



Therefore we assume that these corrections 



can be neglected. Now we can finally discuss higher-order 
corrections to r + (Q 2 ), which represents the ratio of the 
pure "plus" components. Accordingly, we can intepret 
the result in Eq.(5) of Ref. [7[ as higher-order corrections 
to Eq. ((24|) . This interpretation is explicitly confirmed 
up to order a s in Chapter 7 of Ref. [2J, where also the 
same set of equations used in the computation of Ref. [7[ 
are obtained. Further arguments to support it and its 
scheme dependence will be discussed in Ref. [f~6]. We 
denote the approximation in which Eqs. (ll9l27p are used 
as LO + NNLL and the one in which r + (Q 2 ) in Eq. (j2"7) 
is replaced by the result of Eq.(5) in Ref. [7| up to order 
c^J 2 as N 3 LO a p prox + NNLL. That this last one is actu- 
ally a good approximation will be shown below. In both 
approximations considered we can summarize the main 
theoretical result of this Letter in the following way 
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FIG. 1: Gluon and quark multiplicities fits compared to the data. 
The gray dashed line is the LO + NNLL result, the orange solid line 



is the N 3 LO a 



- NNLL result and the red dotted line is the fit 



with four constant coefficients. The orange band corresponds to the 



estimated error of the fitted parameters in the N 3 LO a 



-NNLL 




FIG. 2: Glun-quark multiplicity ratio prediction compared to data. 
The gray solid upper line is the prediction of Ref. Q, the others 
are as in Fig[T] 
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for the gluon-quark multiplicity ratio. Eqs. (|28l29l) de- 
pend only on two parameters, D g (0, Qq) and D s (0, Qq), 
with a simple physical interpretation: they are just the 
gluon and the quark multiplicities at the arbitrary scale 

Qo- 

We have performed a global fit of our resummed formu- 
lae, Eqs. (|28l29p . to the experimental data to extract the 
values of D g (0,Ql) and D s (0,Q 2 ). With Q = 50GeV, 
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90% C.L. 
90% C.L., 

90% C.L. 
90% C.L., 



(31) 



(32) 



the result of the fit is given by 

D g (0,Ql) = 24.31 ±0.85; 
D a (0,Q$) = 15.49 ±0.90; 

in the LO + NNLL case and by 

D g (0,Ql) = 24.02 ±0.36; 
D s {0,Q 2 ) = 15.83 ±0.37; 

in the N 3 LO apP rox + NNLL case in agreement with the 
experimental values within the errors. However, the 90% 
C.L. error in the N 3 LO approx +NNLL case is much smaller 
reflecting a much better fit to the data at all energies. 
Indeed, per degree of freedom we obtain \ 2 — 18.09 in 
the LO + NNLL case, while we have \ 2 — 3.71 in the 
N 3 LO approx + NNLL case. In our analysis, we have used 
the NLO solution for the running coupling according to 
Ea. ([15]) with a s {M z ) = 0.118 and n f = 5. We have 
checked that varying the arbitrary scale Qq does not 
change the resulting value of \ 2 as expected and that 
moving from LL to NNLL the renormalization scale de- 
pendence is strongly reduced. 

In FigfT] we plot the gluon and quark multiplicities 
according to Eqs. (|26l28l29p using the fitted parameters 
given in Eqs. (|31l32[) . Using the data selection of Ref. [2l| , 



the measurements are taken from Refs. [2lM24| for the 
gluon multiplicity and from Refs. 2^, 2(| and references 
therein for the quark multiplicity. The result of a fit 
where the normalization coefficients arc assumed con- 
stant without any additional constraint is also plotted 
showing that N 3 LO ap prox + NNLL is indeed a good ap- 
proximation. To check the consistency of the data sets, 
we have used Eq. (f30|) together with the result of the fit 
from the gluon and quark multiplicities in Eqs. (|31l32[) to 
predict the gluon-quark multiplicity ratio. The result to- 
gether with the corresponding data are shown in Figj2] 
The data are taken from Refs. [H], [12, 25-3(| and refer- 
ences therein, covering essentially all available measure- 
ments. One can see that the data do not agree very well 
at small scales, an isssue that will be discussed elsewhere 

As concluding remarks we remind here that the main 
problem in describing the data was that the theory failed 
badly in the description of the data for the gluon and 
the quark jets simultaneously (or equivalently for the 
ratio r) even if the perturbative series seems to con- 
verge very well. We have shown in this Letter that our 
N 3 LO approx + NNLL result solves this problem explaining 
the discrepancy of the results with the data obtained in 
Ref. 0] as due to the absence of the singlet "minus" com- 
ponent governed by f r _ es {Q, Q 2 , Q§) in Eqs. (j29l30j) . This 
component is included here for the first time. The most 
natural possible future improvement consists in including 
corrections of NLO or beyond to r_(Q 2 ). Our general- 
ized result depends on two parameters, which represent 
our initial condition. They have been fixed performing 



a fit and have a simple physical meaning because they 
just represent the gluon and the quark multiplicity at a 
certain arbitrary scale Qq. We hope that additional mea- 
surements of these observables will come from the LHC 
to test our results on a much wider energy range. 
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